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touchdown.  An  examination  of  three  short  range  prediction  techniques  for 
forecast  times  of  2,  5,  10,  30  and  60  minutes  revealed  that  the  Equivalent 
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SVR  conditions  and  yields  slightly  better  results  than  the  Markov  and  REEP 
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Development  and  Evaluation 
of  a  Tower  Slant 
Visual  Range  System 


1.  INTRODUCTION 

Improved  airfield  low -visibility  capabilities  are  being  achieved  through  the 
installation  of  advanced  instrument  landing  systems  at  military  and  civilian  air¬ 
ports.  With  these  advances,  the  requirement  for  improved  accuracy  in  prediction 
techniques  for  airfield  below-limit  conditions  becomes  more  prominent.  Given  a 
reliable  forecast  technique,  the  operational  capacity  of  an  airfield  can  be  optimized 
and  extended  to  include  Category  Ilia  and  IHb  landings. 

The  Federal  Aviation  Administration  (FAA)1  provides  a  definition  of  Slant 
Visual  Range  (SVR).  Under  daylight  conditions,  SVR  is  the  greater  of  the  slant 
distance  to  (1)  the  farthest  high-intensity  runway  edge  light  or  approach  light  which 
a  pilot  can  see  at  decision  height  (DH)  on  the  approach  path,  or  (2)  the  slant  dis¬ 
tance  which  would  have  a  transmittance  of  5.  5  percent.  Under  conditions  of  dark¬ 
ness,  SVR  is  the  slant  distance  to  the  farthest  high-intensity  runway  edge  light  or 
approach  light  which  a  pilot  will  see  at  DH  on  the  approach  path.  The  requirement 
for  a  glideslope  slant  visual  range  determination  and  short-range  prediction  cap¬ 
ability  exists  within  the  USAF  operating  commands.  Specifically,  this  report  will 


(Received  for  publication  13  September  1979) 

1.  Bradley,  G.  S. ,  Lohkamp,  C.W.,  and  Williams,  R.W.  (1976)  Flight  Test 
Evaluation  of  Slant  Visual  Range/Approach  Light  Contact  Height  (SVR/ 
ALCH)  Measurement  System,  Final  Report  Phase  III,  FAA-RD-76-167. 


will  address  this  scenario  for  Category  I  (DH:200  feet.  Runway  Visual  Range: 

2600  feet)  and  Category  II  (DH:100  feet,  Runway  Visual  Range:  1200  feet)  landing 
operations. 

The  objective  of  this  study  is  to  develop  a  statistically  based  short-range  pre¬ 
diction  technique  which  yields  probability  estimates  of  below -limit  SVR  conditions. 
The  approach  taken  relies  on  the  use  of  atmospheric  extinction  coefficient  meas¬ 
urements  on  a  tower  offset  to  the  side  of  the  landing  zone  similar  in  concept  to  a 

system  proposed  by  the  FAA  and  originally  tested  by  Bradley,  et  al.  1  In  a  pre- 

2 

liminary  analysis,  Hering  and  Geisler  demonstrated  the  importance  of  measure¬ 
ments  above  ground  level  for  a  meaningful  description  of  visibility  conditions  in 
support  of  aircraft  landing  operations.  Measurements  aloft  up  to  DH  were  shown 
to  be  of  special  importance  in  coastal  advection  fog  conditions  since  the  visibility 
along  the  elevated  glide  slope  is  predominately  lower  than  indicated  by  surface  run¬ 
way  visual  range  (RVR)  measurements. 

This  report  documents  the  final  phase  of  the  development  and  evaluation  of  a 
tower  SVR  system  by  AFGL  in  which  the  remote  tower  approach  has  been  extended 
to  include  Category  I  operations.  Three  prediction  techniques  to  yield  probability 
estimates  of  below -limit  SVR  conditions  have  been  evaluated  and  comparisons 
made  to  determine  their  respective  accuracy  and  reliability.  The  first  technique, 
Markov,  is  an  application  of  the  Ornstein-Uhlenbeck  process.  The  second  tech¬ 
nique,  REEP  (Regression  Estimation  of  Event  Probabilities),  is  a  multiple  regres¬ 
sion  analysis  of  a  dichotomous  response  variable  on  dichotomous  independent 
variables.  The  third  technique,  Equivalent  Markov,  is  based  on  a  classical  Markov 
transition  matrix  whose  powers  are  used  to  prepare  forecasts  for  n-time  steps  in 
the  future.  In  the  latter  two  techniques,  probability  estimates  are  direct  outputs 
conditioned  on  the  observed  SVR.  In  the  Markov  technique,  the  observed  SVR 
must  be  transformed  into  an  equivalent  normal  deviate  before  the  probabilities 
can  be  computed.  Forecast  times  examined  in  this  study  are  2,  5,  10,  30  and  60 
minutes. 


2.  TEST  FACILITY  AND  DATA  SETS 

The  development  and  testing  of  a  SVR  system  are  a  part  of  the  field  experi¬ 
ments  at  the  Air  Force  Geophysics  Laboratory  (AFGL)  Weather  Test  Facility  (WTF) 
at  Otis  AFB,  Massachusetts.  The  SVR  technique  development  is  also  part  of  a 


2.  Hering.  W.S.  and  Geisler.  E.B,  (1978)  Forward  Scatter  Meter  Measurement 
of  Slant  Visual  Range,  AFGL-TR-78-0l9l,  AD  A064  425. 
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continuing  program  to  upgrade  the  Modular  Automated  Weather  System  (MAWS) 
developed  by  AFGL  for  fixed  base  Air  Force  requirements. 

Measurements  of  atmospheric  extinction  coefficient  were  obtained  from  for¬ 
ward  scatter  meters  (FSMs)  mounted  on  three  towers  (A.  Q  and  X)  in  the  AFGL 
WTF.  Figure  1  shows  the  configuration  of  the  instrumented  towers  and  the  heights 
at  which  the  FSMs  were  mounted.  The  FSM  instruments  (built  for  the  Air  Force 
by  EG&G)  have  been  used  successfully  in  several  visibility  experiments  carried 
out  by  the  AFGL  over  the  past  several  years.  Analysis  of  the  performance  char¬ 
acteristics  of  the  FSM  for  the  measurement  of  extinction  coefficient  and  sensor 

4  5 

equivalent  visibility  are  given  by  Muench  and  Chisholm  and  Jacobs. 

The  WTF  is  located  in  the  Cape  Cod  area  where  low  visibility  episodes  are 
predominantly  caused  by  heavy  advection  fog  accompanied  by  light  rain  or  drizzle. 
Data  are  collected ’at  a  rate  of  five  observations  per  minute  and  are  placed  on 
magnetic  tape.  Subsequent  processing  yields  a  time  series  of  one-minute  averages 
of  extinction  coefficient  for  each  sensor.  The  continuous  data  stream  from  August 
1977  to  April  1979  was  evaluated  to  identify  reduced  visibility  episodes.  This 
examination  revealed  a  characteristic  difference  between  the  SVR  and  RVR  meas¬ 
urements  during  periods  of  advection  fog  and  rain  of  moderate  or  greater  intensity. 
Figures  2  and  3,  acquired  in  October  1978,  clearly  depict  this.  A  pre-frontal 
band  of  showers  passed  over  the  Cape  Cod  area  and  by  2230  GMT,  heavy  advection 
fog  dominated.  Note  that  during  periods  of  moderate  rain,  the  vertical  gradient  of 
extinction  coefficient  ( 10"*  km"  *)  is  nearly  zero.  During  the  advection  fog  period, 
however,  there  is  a  significant  increase  of  extinction  coefficient  (10_*  km'1)  with 
height.  SVR  specification  equations  developed  from  a  single  data  set  of  these  two 
types  of  low-visibility  episodes  could  reduce  the  probability  of  detecting  a  below- 
limit  condition  of  SVR.  That  is,  a  systematic  bias  towards  the  RVR  measurement 
would  result.  In  order  to  remove  this  systematic  bias,  the  selection  of  low  visibil¬ 
ity  episodes  was  restricted  to  the  advection  fog  type.  Episodes  were  selected  by 
applying  the  criterion  that  the  SVR  in  the  landing  zone  (tower  A)  was  1  km  *  or 
more  for  a  period  of  an  hour  or  greater.  Given  the  universally  accepted  relation¬ 
ships  of  extinction  coefficient  to  visibility,  a  value  of  1  km  1  corresponds  to  1.8  nm 
using  daytime  conditions  (Koschmieder's  Law)  and  3.  8  nm  at  night  (Allard's  Law). 


3.  Tahnk,  W.  R.  and  Lynch,  R.H.  (1978)  The  Development  of  a  Fixed  Base 

Automated  Weather  Sensing  and  Display  System!  AFGL-TR-78-0009. 

AD  A054  805. 

4.  Muench,  H.S.,  Moroz,  E.Y.,  and  Jacobs,  L.  P.  (1974)  Development  and 

Calibration  of  the  Forward  Scatter  Visibility  Meter,  AFCRL-TR-t4-0145, 

AD  785  270. -  - - 


5.  Chisholm,  D.A.  and  Jacobs,  L.  P.  (1975)  An  Evaluation  of  Scattering-Type 
Visibility  Instruments,  AFCRL-TR -75-044 1,  AD  A0l6  766. 


Figure  3.  One  Minute  Average  Values  of  Extinction  Coefficient  as 
Measured  at  Tower  A  for  Advection  Fog  on  14  October  1978 


SVR  specification  and  prediction  equations  were  derived  from  a  dependent  data 
set  which  consisted  of  twenty  reduced  visibility  episodes  and  which  contained  in 
excess  of  6300  minutes  of  data.  The  independent  data  set,  used  to  test  and  compare 
the  various  techniques,  was  drawn  from  twenty-five  other  episodes  totalling  over 
9900  minutes  of  reduced  visibility. 


3.  SLANT  VISUAL  RANGE  SPECIFICATION 
3. 1  Category  II 

2 

In  the  report  on  the  first  phase  of  this  study  (Hering  and  Geisler  ),  the  accu¬ 
racy  of  FSM  measurements  as  specifiers  of  SVR  under  Category  II  conditions  was 
demonstrated.  Those  results  also  indicated  that  during  advection  fog  conditions 
observed  at  Cape  Cod,  a  50  ft  remote  tower  system  yields  SVR  estimates  nearly 
as  accurate  as  a  100  ft  remote  tower  system  in  which  point  visibility  measurements 
are  made  at  the  same  height  as  the  Category  II  DH. 
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The  earlier  study,  applicable  to  Category  11  conditions,  sought  to  compare 
the  ability  of  a  method  which  used  remote  tower  measurements  of  visibility  to 
specify  approach  zone  SVR  (ASVR100)  to  a  method  which  relies  on  the  touchdown 
RVR  measurement.  Table  1  summarizes  the  relevant  statistics  for  the  tower 
option  which  relies  on  measurements  at  50  ft  (referred  to  as  Method  2  in  reference 
2),  identified  here  as  ASVR100,  and  the  control  technique  (RVR).  The  threshold 
5  km  *  corresponds  to  about  1  2  mile  (800  m)  daytime  visibility  and  1  mile 
(1600  m)  at  night.  The  threshold  12  km  *  corresponds  to  1/4  mile  (400  m)  daytime 
and  1/2  mile  (800  m)  at  night.  Appendix  A  describes  the  method  of  calculating  the 
statistics  used  here.  Clearly,  aSvr1^30  yields  superior  specifications  of  ASVR100 
than  does  RVR.  Figure  4,  which  is  a  time  series  plot  of  an  advection  fog  episode 
in  October  1978,  demonstrates  that  the  major  deficiency  of  the  RVR  method  lies  in 
the  persistent  optimism  it  conveys  to  the  pilot  at  DH  of  seeing  his  reference  point. 
Recall  that  lower  extinction  coefficient  values  reflect  higher  visibility. 

3.2  Category  I 

The  measurement  of  SVR  to  200  feet  was  made  possible  with  the  installation  of 
FSMs  on  tower  A  at  the  150  and  200  ft  levels.  The  objective  here  is  to  develop  a 
specification  equation  for  SVR  of  200  feet  (ASVR200)  in  the  simulated  approach 
zone  using  FSM  measurements  of  extinction  coefficient  taken  at  towers  Q  and  X. 

The  FSM  measurements  at  discrete  points  in  the  vertical  are  converted  to  a 
weighted  vertical  average  which  is  used  to  represent  SVR  in  the  approach  zone 
through: 

ASVR200  -  (A200  +  2A150  +  2A100  *•  2A50  ♦  A10)  8  .  (1) 

Application  of  the  multiple  linear  regression  technique  to  the  dependent  data 
set  of  6336  observations  of  reduced  visibility  yielded  two  equations  for  specifying 
approach  zone  SVR  (ASVR200).  The  first  equation  specifies  ASYR200  from  remote 
tower  (Q10Q  and  Q 50)  and  touchdown  RVR  (X10)  measurements  of  extinction  coef¬ 
ficient  as  follows: 

EST1  -  Q100  +  0.  86X10  -  0.  61Q50  (2) 

The  second  equation  specifies  ASVR200  from  remote  tower  (Q50)  and  touch¬ 
down  RVR  (X10)  measurements  of  extinction  coefficient  as  follow  s: 

EST2  1.06Q50  •  0.547X10  (3) 

The  second  equation  (EST2)  was  developed  to  measure  the  loss  in  skill  (if  any)  which 
would  be  experienced  if  only  a  50  ft.  offset  tower  could  be  utilized. 
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Table  1.  Statistical  Results  of  ASVR100  and  RVE  as  Specifiers 
for  ASVH100 


Statistic 

aSvrToo 

RVR 

Correlation  Coefficient 

0.  979 

0.  966 

Per  cent  RMSE 

25.  4 

40.  0 

Rias 

0.4 

-37.  9 

Threshold:  5  km  * 

TS 

90 

71 

POD 

93 

71 

FAR 
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Figure  4.  One  Minute  Average  Values  or  ASVR100,  ASVR100  and  RVR  for 
Advert  ion  Fog  on  14  October  1978 
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Table  2.  Statistical  Results  of  EST1,  EST2  and  RVE  as 
Specifiers  for  ASVR200 


Statistic 

EST1 

EST2 

RVR 

Correlation  Coefficient 

0.  971 

0.  931 

0.899 

Per  cent  RMSE 

30.  8 

44.  1 

58.  9 

Bias 

-5.  6 

-11.  2 

-82.9 

Threshold:  5  km"' 

TS 

89 

68 

55 

POD 

91 

69 

55 

FAR 

2 

1 

0 

Threshold:  12  km'1 

TS 

87 

70 

32 

POD 

91 

71 

32 

FAR 

5 

1 

0 

Table  2  summarizes  the  relevant  statistics  for  these  specification  equations 
and  the  control  technique  (touchdown  RVR). 

Unlike  the  results  under  Category  II  conditions,  the  statistics  of  Category  I 
conditions  show  that  the  100  ft  method  (EST1)  outperforms  the  50  ft  method 
(EST2)  by  nearly  15  percentage  points  in  the  RMSE.  Similar  differences  exist  in 
the  other  statistics.  Again  we  find  the  remote  tower  equations  are  superior  to  the 
touchdown  RVR  measurement.  To  demonstrate  these  results  we  select  ASVR200, 
RVR  and  ASVR200  (redesignated  EST1)  and  display  them  in  a  time  series  of  an 
advectionjog  episode  in  Figure  5.  Note  the  stubborn  bias  of  the  RVR  measurement 
while  ASVR200  closely  tracts  ASVR200. 


4.  PREDICTION  TECHNIQUES 


4.1  Markov  Model 

g 

Gringorten0  adapted  for  meteorological  use  a  special  class  of  the  Markov  chain 
called  the  Ornstein-Uhlenbeck  process  in  which  a  single  continuous  variate  serves 

7 

as  both  predictor  and  predictand  (through  a  time  lag).  Bering  and  Quick  used  this 

model  with  good  success  in  forecasting  atmospheric  extinction  coefficient  in  the 

AFGL  Mesonetwork  experiments.  This  model  was  also  used  as  a  control  technique 

in  the  evaluation  of  more  complex  forecasting  techniques  appropriate  for  radiation 
g  a 

fog  (Tahnk  )  and  advective  situations  (Chisholm  ). 

In  the  Ornstein-Uhlenback  stochastic  process,  the  value  St  of  the  continuous 

variate  at  time  t  is  related  to  an  initial  value  S  as  follows: 

o 


St  =  V‘o+  t1  -<%)2' 


(4) 


where  a  =  one  hour  autocorrelation  coefficient:  t  =  forecast  time  interval  (hours); 
o 

p  -  normalized  probability. 


6.  Gringorten,  1. 1.  (1972)  Conditional  probability  for  an  exact  noncategorized 

initial  condition,  Mon,  Wea.  Rev.  100:796-798. 

7.  Hering,  W.S.  and  Quick,  D.  L.  (1974)  Hanscom  visibility  forecasting  exper¬ 

iments,  Proc,  of  5th  Conference  in  Weather  Forecasting  and  Analysis, 
American  Meteorological  Society,  pp.  224-227, 

8.  Tahnk,  W.R.  (1975)  Objective  Prediction  of  Fine  Scale  Variations  in  Radiation 

Fog  Intensity,  AFCRL-TR -75-0269,  AD  A0l4  774. 

9.  Chisholm,  D.A.  (1976)  Objective  Prediction  of  Mesoscale  Variations  of  Sen¬ 

sor  Equivalent  Visibility  During  Advection  Situations.  AFGL-TR-76-0132. 
AD  A030  332. 
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During  the  AFGL  Mesonetwork  experiment,  pQ  was  evaluated  in  an  effort  to 
increase  the  accuracy  of  the  model  when  using  high-frequency  (one-minute)  obser¬ 
vations.  The  evaluation  yielded  seasonal  values  of  pQ  =  0.96  for  the  period 
November  to  March  and  p  =  0.  93  from  April  to  October. 

SQ  is  determined  in  the  following  way: 

SQ  =  k(lnSVR0)  +1  (5) 

where 

InSVR  =  Slant  Visual  Range  (Extinction  Coefficient)  at  time  t  in  logarithmic 
0  form  ° 

k,  1  =  function  of  time  before  or  after  sunrise  (Ats)  through 

k  =  aAts  +  b  (6) 

1  =  cAts  +  d  (7) 

where 

a,  b,  c,  d  =  constants  and  coefficients 

Ats  =  hours  away  from  sunrise  . 

The  unconditional  cumulative  frequency  distribution  (cfd)  of  prevailing  visibil¬ 
ity  defined  by  the  terms  a,  b,  c  and  d  is  dependent  on  season  and  time  of  day.  They 
were  determined  from  twenty-three  years  of  hourly  observations  made  by  opera¬ 
tional  observers  at  Otis  AFB's  runway  adjacent  to  the  AFGL  WTF.  An  assumption 
is  made  which  equates  the  cfd  of  prevailing  visibility  to  the  cfd  of  SVR.  This  is  a 
necessary  assumption  because  there  does  not  exist  a  SVR  data  base  to  accurately 
determine  its  cfd. 

SVRq  for  Categories  I  and  II  conditions  were  obtained  from  the  regression 
equations  which  yield  specifications  of  ASVR200  and  ASVR100  respectively.  Once 
determining  St  from  Sq  and  the  appropriate  pQ,  one  computes  the  probability  of 
exceeding  certain  operationally  significant  thresholds  by  assuming  a  Gaussian  sto¬ 
chastic  process  and  using  an  approximation  to  the  normal  probability  integral. 

In  summary,  the  most  probable  normalized  SVR  is  first  determined  based  on 
a  Markov  decay  of  the  initial  SVR.  Assuming  a  Gaussian  distribution  of  errors 
about  the  most  probable  SVR,  one  is  able  to  calculate  the  probability  of  exceeding 
a  threshold. 
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4.2  Regreaion  Estimation  of  Event  Probabilitie*  (REEF) 

REEP  calculates  probabilities  of  being  within  categories  which  can  easily  be 
converted  into  exceedance  probabilities.  Instead  of  transforming  the  initial  SVR 
into  a  most  probable  SVR  from  which  an  exceedance  probability  can  be  found 
(Markov  Model),  REEP  uses  the  initial  SVR  directly. 

REEP10  is  similar  to  a  forward  screening  regression  technique  where  a  sub¬ 
set  of  predictors  is  selected  from  a  large  set  of  possible  predictors.  The  relation¬ 
ship  between  the  predictand  and  the  selected  predictors  is  a  linear  function  whose 
coefficients  are  determined  by  least  squares. 

REEP  requires  that  each  raw  variable  be  broken  down  into  a  set  of  two  or 
more  mutually  exclusive  and  exhaustive  categories  called  dummy  variables.  The 
dummy  variable  categories  are  assigned  a  value  of  1  if  the  continuous  variable 
falls  within  the  range  of  the  category  limits,  or  0  if  it  is  outside  the  range.  There¬ 
fore,  for  any  one  observation,  one  dummy  variable  will  be  assigned  a  value  of  l 
and  all  of  the  others  will  be  assigned  a  value  of  0. 

The  result  of  the  least  squares  technique  is  a  set  of  equations  one  for  each 
predictand  category: 


P 

P 


1 

2 


=  A01  +  AUX1 
=  A02  +  A12X1 


+  A21X2  + 


+  A22X2  + 


+  A91X9 
+  A92X9 


(8) 


P5  =  A05  +  A 15X1  +  A25X2  +  ’  ’ *  +  A99X9 

Five  predictand  category  limits  were  selected  to  coincide  with  the  thresholds 
prescribed  beforehand  and  which  are  listed  in  Table  3. 

Category  limits  were  assigned  to  each  predictor  based  on  the  relative  fre¬ 
quency  distribution  of  the  SVR  values  in  the  dependent  data  set.  ASVR200 
was  used  to  establish  the  frequency  distribution  of  the  dummy  predictors  for 
Category  I  conditions.  Similarly,  ASVR100  was  used  for  Category  II  conditions. 
The  resultant  categories  for  the  predictors  are  listed  in  Table  4. 

REEP  is  formulated  to  insure  internal  consistency  among  predictand  cate¬ 
gories  such  that  the  sum  of  probabilities  totals  unity.  This  is  insured  by  using 


10.  Miller,  R.G.  (1964)  Regression  Estimation  of  Event  Probabilities,  Tech.  Rpt. 
7411-121,  The  Travelers  Research  Center,  Inc.,  Hartford,  Conn. 
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Table  3.  REEP  Predictand  Category  Limits 


Predictand  SVR  Limits 
(km”1) 

Predictand  Category 

Greater  Than 
or  Equal  To 

Less  Than 

1 

1.0 

2.  0 

2 

2.0 

5.0 

3 

5.  0 

12.  5 

4 

12.  5 

28.  1 

5 

28.  1 

80.0 

Table  4.  REEP  Predictor  Category  Limits 


Predictor  SVR  Limits 
(km”1) 

Predictor  Category 

Greater  Than 
or  Equal  To 

Less  Than 

1 

1.0 

2.0 

2 

2.0 

3.  5 

3 

3.  5 

5.0 

4 

5.0 

8.  5 

5 

8.  5 

12.  5 

6 

12.  5 

18.  5 

7 

18.  5 

25.0 

8 

25.  0 

40.0 

9 

40.  0 

80.  0 

the  same  predictor  categories  for  each  predictand  equation.  Each  predictand 
equation  yields  a  probability  of  occurrence  for  the  category  (P.).  REEP  prediction 
equations  (for  intervals  of  2,  5,  10,  30  and  GO  minutes)  were  generated  from  a 
randomly  selected  sub-sample  of  3000  observations  from  the  dependent  sample. 
The  order  of  selection  of  dummy  variable  predictors  and  the  REEP  equation  sets 
are  listed  in  Appendix  B. 

Conversion  of  REEP  categorical  probabilities  to  the  probability  of  exceeding  a 
threshold  is  simply  made  after  the  P.,  i  =  1,  2,  ...  5,  are  calculated.  For 
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example  the  probability  of  SVR  being  below  1  !2  mile  daytime  (P^  or  greater)  given 
Pj  =  0.  1,  P2  0.  15,  P3  -  0.35,  P4  -  0.  25  and  Pg  0.  15  is 

P(SVR  <1/2  mile  day)  P(SVR  >  5  km'1)  =  P3  i  P4  f  P. 

Hence 

P(SVR  <  1/2  mile  day)  =  0.  75  . 

4.3  Equivalent  Markov  Model  Baaed  on  RF.F.P 

In  principle  it  is  possible  to  develop  a  prediction  scheme  in  which  SVR  changes 
are  modeled  as  a  discrete  Markov  process  with  a  finite-state  space.  The  classical 
Markov  transition  Matrix  St  can  be  used  in  preparing  forecasts  for  any  number  of 
time  steps  (n)  in  the  future  by  using  powers  of  X!.  A  prediction  method  that  yields 
probabilistic  forecasts  comparable  to  the  classical  Markov  process  but  without 
the  necessity  of  utilizing  S/I  explicitly  has  been  proposed  by  Miller.  Miller  pro¬ 
posed  using  the  coefficients  from  a  set  of  REEP  equations  to  determine  the  one- 
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step  transition  matrix.  Whiton  found  this  technique  to  yield  results  comparable 

to  the  classical  Markov  process  but  much  easier  to  develop  and  apply  to  practical 

13 

forecasting  problems.  Hillier  and  Lieberman  utilize  a  set  of  conditional  or 
transition  probabilities  to  define  the  one-step  transition  matrix  P.  It  is  this  latter 
approach  that  the  technique  used  here  is  based  upon. 

The  REEP  equations  provide  a  conditional  probability  matrix  P  equivalent  to 
M.  One  assumes  that  the  categories  (predictor  and  predictand)  used  in  REEP  form 
a  finite  Markov  chain  (Xj),  i  =  1,  2,  ....  5.  Hence,  P  is  a  square  matrix  whose 
predictor  and  predictand  category  limits  are  those  found  in  Table  3. 

The  stochastic  process  (X^)  has  the  Markovian  Property 

P(Xt+1  i  !  Xt  --  ;i)  --  P..  :  i,j  =  1,  2 . 5  .  (9) 

The  existence  of  one-step  transition  probabilities  P^  implies  that  for  each  i,  j 
and  n 


11.  Miller,  R.G.  (1968)  A  Stochastic  Model  for  Real-Time  On-Demand  Weather 

Predictions,  Proc.  1st  Statistical  Meteorological  Conference,  American 
Meteorological  Society,  pp.  48-51, 

12.  Whiton,  R.W.  (1977)  Selected  Topics  in  Stat istical  Meteorology,  AWS-TR- 

77-273,  pp.  7-1  to  7-45  or  Chap.  7  Markov  Processes.  ** 

13.  Hillier,  F.S.  and  Lieberman,  G.J.  (1974)  Operations  Research,  Holden-Day, 

Inc.,  pp.  351-369. 
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P(Xt+n  =  i  I  Xt  -  j)  =  HXn  =  i  I  XD  -  j)  =  P4j 


(n) 


<10) 


P.  ^  is  the  conditional  probability  that  the  random  variable  X,  starting  in  state  i, 
will  be  in  state  j  after  exactly  n  steps  (time  units)  and  is  called  the  n-step  transi¬ 
tion  probability. 

Since  the  p.^ are  conditional  probabilities  they  must  have  the  properties 


p.  >  0 


for  all  i,  j  =  1,  2,  ....  5 


(11) 


ZP.  .(n>  =  1  for  all  i  =  1,  2 . 5 

ij 

3  =  1 


(12) 


is  defined  as  the  n-step  transition  matrix  whose  elements  are  P. 


p  (n)  „  (n) 

*11  *12 


(n) 
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p  (n)  „  (n)  p  (n) 

”51  *52  *  •  *55 


(13) 


The  Chapman -Kolmogorov  equations  provide  a  method  for  computing  these 
n-step  transition  probabilities: 


pi3<n>  =  I  Pik<V>  Pk/n"V>  -  for  all  i.j,nandl<v<n 


(14) 


k=l 


More  generally,  it  follows  that  the  matrix  of  n-step  transition  probabilities 
can  be  obtained  from  the  expression: 

P(n)  --  P.P...P  -  P"  =  PPn_1  -  Pn_1  P  .  (15) 

Thus  the  n-step  transition  probability  matrix  can  be  determined  by  computing 
the  nth  power  of  one-step  transition  matrix  P. 

In  order  to  calculate  the  probability  of  exceeding  a  threshold  one  needs  to  first 
pre -multiply  P°  by  the  observation  matrix  (f: 
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(16) 


*  =  (011012013014°15> 


where  the  0.^  category  is  identical  to  those  listed  in  Table  3.  Having  calculated 
5Pn  one  computes  the  exceedance  probabilities  in  the  same  way  as  in  the  REEP 
technique. 

In  the  initial  application  of  the  Equivalent  Markov  model,  a  two  minute  one- 
step  transition  matrix  was  utilized  to  generate  predictions  for  b  through  60  min- 
uites  and  comparisons  were  made  to  Markov  and  REEP  results.  The  Equivalent 
Markov  model  performed  poorly  when  applied  in  this  way,  especially  at  the  longer 
forecast  intervals.  In  an  attempt  to  rectify  this,  a  variation  of  the  Equivalent 
Markov  technique  was  developed  and  tested  wherein  consideration  was  given  to  the 
use  of  two  separate  transition  matrices  to  cover  the  full  range  of  predictions. 

For  SVR  of  100  feet  (Category  n  operations),  the  best  results  were  obtained  with 
a  two  minute  one-step  transition  matrix  to  compute  conditional  probabilities  for 
forecasts  of  five  minutes  or  greater.  For  SVR  of  200  feet  (Category  I  operations), 
a  two  minute  one-step  transition  matrix  is  used  to  compute  two  and  five  minute 
forecasts  and  a  ten  minute  one-step  transition  matrix  is  used  to  compute  forecasts 
of  ten  minutes  or  greater.  (The  calculation  of  the  five  minute  forecast  is  made  by 
calculating  a  simple  average  of  the  four  and  six  minute  forecasts. ) 


Table  5.  Values  of  Percent  Improvement  of  P-Score  for  Equivalent 
Markov  and  Variation 


SVR  Height 
(ft) 

Forecast  Time 
(min) 

Thresholds  (mile) 

1 

1/2 

1/4 

200 

5 

65 

35 

49 

10 

37 

37 

36 

30 

43 

34 

-2 

60 

33 

26 

-18 

100 

5 

43 

30 

62 

10 

49 

39 

62 

30 

54 

42 

53 

60 

49 

38 

44 

Table  5  presents  the  percent  improvement  in  P-score  (defined  in  Section  5)  of 
the  refined  approach  over  the  original  approach  when  evaluated  on  the  dependent 
data  set.  Clearly  using  two  one-step  transition  matrices  is  a  substantial 
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Improvement  over  using  one  matrix  and  this  variation  alone  will  be  applied  In  the 
rest  of  the  study.  The  transition  matrices  can  be  found  In  Appendix  C. 

5.  EVALUATION  CRITERIA 

It  Is  desirable  that  statements  of  the  probability  of  a  weather  event  by  reliable; 
that  is.  over  a  period  of  time  the  event  should  actually  occur  with  the  frequency 
Implied  by  the  probability  forecast.  It  Is  also  desirable  that  there  be  resolution  in 
the  probabilities;  that  is,  they  should  be  as  close  to  0  or  to  100  as  possible  when 
the  event  does  not  occur  or  does  occur,  respectively. 

The  P-score  as  defined  by  Epstein**  and  Murphy*^  measures  these  two  char¬ 
acteristics  of  probability  forecasts  and  represents  the  mean  square  difference 
between  the  forecast  and  observed  probability  distribution.  It  is  given  by: 

N 

PS  =  ^  I  <Ft  -  0.)2  (17) 

i=  1 

where 

Fj  =  forecast  probability  of  exceeding  the  threshold 

Ch  =  verification  based  on  observed  SVR  as  follows: 

1  =  if  threshold  is  exceeded 
0  =  if  not  . 

The  nature  of  PS  is  such  that  it  attains  a  value  of  0  with  perfect  probability 
forecast  (that  is,  0  or  100  percent  probability  and  correct)  or  a  value  of  1  if  totally 
imperfect  (that  is,  0  or  100  percent  and  incorrect). 

6.  RESULTS:  PREDICTION  TECHNIQUES 

The  aviation  industry  has  progressed  to  the  point  where  a  significant  number 
of  airfields  have  been  upgraded  and  approved  for  Category  II  operations.  As  such, 
this  study  has  focussed  on  Categories  I  and  II  landing  conditions.  These  categories 

14.  Epstein,  E.  S.  (1969)  A  scoring  system  for  probability  forecasts  of  ranked 

categories,  J.  Appl.  Meteorol.  8:985-987. 

15.  Murphy,  A.H.  (1969)  On  the  ranked  probability  score,  J.  Appl.  Meteorol. 

8:988-989.  - 
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correspond  to  a  daytime  SVR  threshold  of  1/2  mile  for  Category  I  and  a  threshold 
1/4  mile  for  Category  II,  respectively.  The  three  prediction  techniques  were 
applied  to  each  of  the  twenty-five  independent  data  episodes.  The  results  of  each 
episode  were  combined  to  reflect  the  overall  accuracy  and  reliability  of  each  pre¬ 
diction  technique.  Comparison  of  the  three  techniques  was  achieved  by  use  of  the 
P-score  statistic  which  calculates  the  mean  squared  probability  error  and  by 
reliability  graphs. 

Verification  results  of  the  P-score  for  Category  I  operations  given  in  Table  6 
reflect  that  the  Equivalent  Markov  technique  yields  a  slight  improvement  over  the 
Markov  and  REEP  techniques  (2  percent  and  9  percent  respectively).  Figures  6 
to  8  depict  the  reliability  of  the  three  techniques  for  10  min  predictions.  Clearly 
the  predicted  probabilities  of  the  Equivalent  Markov  and  REEP  techniques  are  quite 
close  to  the  observed  frequencies.  However,  the  predicted  probabilities  of  the 
Markov  technique  tend  to  underestimate  the  observed  frequencies. 

Figures  9  and  10  reveal  a  slight  tendency  of  the  Equivalent  Markov  technique 
to  overestimate  below  limit  SVR  conditions  for  Category  I  operations  for  forecast 
times  of  30  to  60  minutes,  respectively. 

Table  7  summarizes  the  results  for  Category  II  operations.  The  Equivalent 
Markov  technique  yields  systematically  better  P-scores  over  the  Markov  and 
REEP  techniques.  Figures  11  and  13  show  the  excellent  correspondence  between 
the  predicted  probabilities  and  the  observed  frequencies  for  the  Equivalent  Markov 
and  REEP  techniques.  Meanwhile,  Figure  12  reflects  the  tendency  of  the  Markov 
technique  to  overestimate  the  observed  frequencies.  A  dropoff  in  skill  of  the 
Equivalent  Markov  technique,  as  the  forecast  interval  is  increased  from  10  to 
60  minutes,  is  due  in  part  to  decreasing  resolution  in  the  low  probability  forecasts. 
This  can  be  seen  in  Figures  11,  14  and  15  that  show  the  percentage  of  forecasts 
in  the  0-5  percent  range  lowers  from  85  to  72  to  58  percent. 


Table  6.  Values  of  P-Score  Category  I  Operations 
Equivalent  Markov,  Markov  and  REEP  Techniques 


Prediction  Technique 

Lag 

Markov 

REEP 

Equiv-Markov 

2 

0.046 

0.  048 

0.  045 

5 

0.  055 

0.  060 

0.  03  5 

10 

0.075 

0.081 

0.  074 

30 

0.  131 

0.  139 

0.  129 

60 

0.  172 

0.  196 

0.  169 
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PREDICTED  PROBABILITY  OF  PREDICTED  PROBABILITY  OF 


CATEGORY  I 
EQUI V-ttARKQV 
ID  MIN  FORECAST 


OBSERVED  RELATIVE  NO.  OF 

FORECASTS!  7. ) 

Figure  6.  Reliability  Graph  for  Category  I  Operations: 

Equivalent  Markov  Technique:  10  Minute  Forecast 


CATEGORY  I 
MARKOV 

10  MIN  FORECAST 


OBSERVED  RELATIVE  NO.  OF 

FORECASTS!  7.) 

Figure  7.  Reliability  Graph  for  Category  I  Operations: 

Markov  Technique:  10  Minute  Forecast 
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PREDICTED  PROBABILITY  OF  PREDICTED  PROBABILITY  OF 

BELOW-LIMIT  VISIBILITY  BELOW-LIMIT  VISIBILITY 


CATEGORY  I 
REEP 

10  MIN  FORECAST 


OBSERVED  RELATIVE  NO.  OF 

FORECASTS!  7. ) 

Figure  8.  Reliability  Graph  for  Category  I  Operations: 

REEP  Technique:  10  Minute  Forecast 

CATEGORY  I 
EQUI V-MARKOV 
30  MIN  FORECAST 


OBSERVED  RELATIVE  NO.  OF 

FORECASTS!  ’/. ) 

Figure  9.  Reliability  Graph  for  Category  I  Operations: 

Equivalent  Markov  Technique:  30  Minute  Forecast 
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PREDICTED  PROBABILITY  OF 
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60  MIN  FORECAST 


20  40  60  80  100 
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RELATIVE  NO.  OF 
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Figure  10.  Reliability  Graph  for  Category  I  Operations: 
Equivalent  Markov  Technique:  60  Minute  Forecast 

CATEGORY  II 
EQU I V-MARKOV 
10  MIN  FORECAST 


OBSERVEO 


RELATIVE  NO.  OF 
FORECASTS! '/.  1 


Figure  11.  Reliability  Graph  for  Category  II  Operations: 
Equivalent  Markov  Technique:  10  Minute  Forecast 


28 


PREDICTED  PROBABILITY  OF  PREDICTED  PROBABILITY  OF 


CATEGORY  II 
EQUIV-MARKOV 
30  MIN  FORECAST 


RELATIVE  NO.  OF 
FORECASTS!  7. ) 


Figure  14.  Reliability  Graph  for  Category  II  Operations: 
Equivalent  Markov  Technique:  30  Minute  Forecast 


CATEGORY  II 
EQUIV-MARKOV 
60  MIN  FORECAST 


OBSERVED  RELATIVE  NO.  OF 

FORECASTS! X ) 

Figure  15.  Reliability  Graph  for  Category  n  Operations; 
Equivalent  Markov  Technique:  60  Minute  Forecast 
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Table  7.  Values  of  P-Score  Category  II  Operations 
Equivalent  Markov,  Markov  and  REEP  Techniques 


Prediction  Technique 

Lag 

Markov 

REEP 

Equiv-Markov 

2 

BSI 

wmmmm 

5 

BEllnB 

WSSM 

10 

30 

0.052 

n 

■ 

60 

0.067 

0.056 

|| 

It  was  stated  above  that  the  Markov  technique  underestimates  the  observed 
frequencies  for  Category  I  conditions  and  overestimates  them  for  Category  II 
conditions.  This  probably  reflects  a  weakness  in  the  original  assumption  of  equat¬ 
ing  the  cfd  of  prevailing  visibility  to  the  cfd  of  SVR,  an  assumption  which  had  to  be 
made  due  to  the  lack  of  SVR  data  to  determine  its  cfd  directly. 


7.  CONCLUSIONS 


Analysis  of  the  data  collected  at  the  AFGL  WTF  demonstrates  the  accuracy 
and  reliability  of  a  remote  tower  SVR  system.  For  both  Categories  I  and  II 
operations  such  a  system  would  provide  a  probability  of  detection  of  below -limit 
conditions  of  91  percent  or  greater,  a  large  improvement  over  using  the  surface 
RVR  instrument.  In  an  examination  of  the  short  range  prediction  techniques,  the 
Equivalent  Markov  technique  was  shown  to  provide  accurate  and  reliable  forecasts 
of  below -limit  SVR  conditions  and  yield  slightly  better  results  on  independent  data 
than  did  the  Markov  and  REEP  technique.  Due  to  the  high-frequency  fluctuations 
characteristic  of  coastal  advection  fog  episodes,  a  continuous  updating  and  dissem¬ 
inating  system  is  recommended  to  allow  this  information  to  be  used  in  a  timely 
manner.  An  automated  weather  system  such  as  MAWS  can  easily  accommodate 

this  feature.  _ 

The  SVR  specification  equations  (XsVRIOO  and  ASVR200)  and  the  Equivalent 
Markov  technique  are  being  integrated  into  the  MAWS  recently  installed  at  the 
AFGL  WTF  as  part  of  a  continuing  investigation  of  automated  airfield  weather 
observation  and  prediction  systems. 
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Fi  =  ASVR200 
Ol  =  ASVR200 
BLAS  =  0  is  perfect 
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3.  Correlation  Coefficient  (R) 


r  -  SXY 

v/fEX^HEY5? 


where 


X  =  (ASVR200  -  ASVR200 


Y  =  (ASVR200  -  ASVR200) 


4.  Contingency  Table 

The  following  data  display  in  Table  Al  depicts  the  accuracy  of  ASVR200  as  a 
specifier  for  ASVR200  for  a  threshold  of  5  km*1  (day:  1/2  mile:  night:  1  2  mile) 
and  a  time  lag  of  0  minutes. 

1)  Threat  Score  (TS) 

TS  =  (5453/(6028  +  97))  X  100  =  89.0% 

2)  Probability  of  Detection  (POD) 

POD  =  (5453/6028)  X  100  =  90.5% 

3)  False  Alarm  Ratio  (FAR) 

FAR  =  (97/5550)  X  100  =  1.7% 


Table  Al.  Contingency  Table  for  Independent  Data  Set  Which 
Illustrates  Method  Computer  Verification  Scores 


Forecast 

Observed 

SVR  =»  5  km*1 

SVR  <  5  km"1 

Total  Min 

Visibility 
below  Limit 

SVR  >5  km-1 

5453 

575 

6028 

Visibility 
above  Limit  . 
SVR  <  5  km"1 

97 

3782 

3879 

Total  Min 

5550 

4357 

9907 
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Appendix  B 

REEP  Technique  Equations 


Category  I  Operations 


Two  Minute  Forecast 

Order  (i) 

P1 

P2 

P3 

P4 

P5 

Predictor 

Xi 

Limits  of 
Xi 

0 

0 

0.  100 

0.025 

0.  667 

0.  198 

0 

1 

0 

-0.  096 

-0.019 

-0.511 

0.627 

8 

25-40 

2 

0 

-0.  100 

-0.025 

-0.671 

0.795 

9 

40-80 

3 

0.  836 

0.  035 

0.005 

-0.677 

-0. 198 

1 

1-2 

4 

0.499 

0.314 

0.058 

-0.672 

-0.  198 

2 

2-3.  5 

5 

0.033 

0.014 

0.  693 

-0. 542 

-0. 198 

4 

5-8.  5 

6 

0.025 

-0.010 

0.465 

-0.281 

-0.  198 

5 

8.  5-12.  5 

7 

0.  122 

0.213 

0.  515 

-0.  651 

-0. 198 

3 

3.  5-5 

8 

0.019 

-0.029 

0.  082 

0.  126 

-0.  198 

6 

12.  5-18.  5 
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Category  I  Operations  (Cont. ) 


Five  Minute  Forecast 

Order 

P1 

P2 

P3 

P4 

P5 

Predictor 

Xi 

Limits  of 
Xi 

0 

0.  008 

0.  102 

0.03 

0.644 

0.  216 

0 

1 

-0.  008 

-0.  098 

-0.  021 

-0.468 

0.  595 

8 

25-40 

2 

-0.008 

-0.  102 

-0.  03 

-0.  631 

0.771 

9 

40-80 

3 

0.  699 

0.  132 

0.  029 

-0.  644 

-0.216 

1 

1-2 

4 

0.482 

0.304 

0.050 

-0. 620 

-0.216 

2 

2-3.  5 

5 

0.016 

0.040 

0.  655 

-0.494 

-0.216 

4 

5-8.  5 

6 

0.017 

-0.  014 

0.419 

-0. 205 

-0.216 

5 

8.  5-12. 5 

7 

0.  122 

0.  194 

0.  510 

-0. 609 

-0.216 

3 

3.  5-5 

8 

0.006 

-0.  011 

0.077 

0.  134 

-0.  205 

6 

00 

H 

1 

m 

r*( 

Order 


Order  (i) 


Ten  Minute  Forecast 


Predictor  Limits  of 
Xi  Xi 


Thirty  Minute  Forecast 


-0.031 
0.  164 
0.043 


2  3 


470 
446 
431 
433 
290 
301 

0.008  -0.295 

0.184  -0.131 

0.019  0.004 


4  *5 


0.364  0.044 

-0.126  0.661 
-0.090  0. 643 

0.  147  0.292 

-0.299  -0.033 

-0.248  -0.044 

0.244  0.074 

-0.208  -0.009 

-0.046  -0.020 


Predictor  Limits  of 
Xi  Xi 


25-40 
40-80 
18.  5-25 
2-3.5 
1-2 

12.  5-18.  5 
3.  5-5 
5-8.  5 


Sixty  Minute  Forecast 


Order  (i) 

P1 

P2 

P3 

P4 

P. 

J 

Predictor 

Xi 

Limits  of 
Xi 

0 

0.066 

0.093 

0.381 

0.339 

0.  120 

0 

1 

-0.048 

-0.060 

-0.343 

0.003 

0.447 

8 

25-40 

2 

-0.009 

-0.024 

-0.327 

0.  141 

0.  219 

7 

18.  5-25 

3 

0.  170 

0.  241 

-0. 113 

-0.  197 

-0.  101 

2 

2-3.  5 

4 

0.  209 

0.  223 

-0.  153 

-0.  211 

-0.067 

1 

1-2 

5 

0.021 

0.  211 

-0.  051 

-0. 157 

-0.024 

3 

3.  5-5 

6 

-0.  007 

-0.034 

-0.  231 

-0.026 

0.298 

9 

40-80 

7 

-0.  047 

-0.  022 

-0.  138 

0.066 

0.  140 

6 

12.  5-18.  5 

8 

-0.  019 

0.  047 

-0.048 

-0.  056 

0.  076 

5 

8.  5-12.  5 

Category  II  Operations 


Two  Minute  Forecast 

Order  (i) 

P1 

P2 

P3 

P4 

P5 

Predictor 

Xi 

Limits  of 
Xi 

0 

0.  869 

0.  114 

0.000 

0 

1 

0.  907 

-0.  865 

-0. 114 

0.000 

1 

1-2 

2 

-0.  869 

0.032 

0.854 

9 

40-80 

3 

-0.  107 

-0.  855 

0.  207 

0.  666 

8 

25-40 

4 

-0.017 

-0.  817 

0.  832 

0.003 

7 

18.  5-25 

5 

0.000 

-0.  736 

0.753 

0.000 

6 

12.  5-18.  5 

6 

0.251 

0.  688 

-0.825 

-0. 114 

0.000 

2 

2-3.  5 

7 

0.  637 

-0.  564 

-0. 114 

0.000 

3 

3.  5-5 

8 

0.019 

0.  223 

-0.  141 

-0.  101 

0.000 

4 

5-8.  5 

Five  Minute  Forecast 

Order  (i) 

P1 

P2 

P3 

P4 

P5 

Predictor 

Xi 

Limits  of 
Xi 

0 

0.  790 

0.  162 

0 

1 

0.  859 

-0.773 

-0.  162 

1 

1-2 

2 

-0.  784 

-0. 021 

9 

40-80 

3 

-0.  776 

0.  183 

8 

25-40 

4 

0.  000 

-0.  048 

-0.  722 

0.  766 

0.005 

7 

18.  5-25 

5 

0.000 

-0.  048 

-0. 623 

0.672 

0.000 

6 

12.  5-18.  5 

6 

0.  280 

0.  583 

-0.704 

-0.  159 

0.000 

2 

2-3.  5 

7 

0.  074 

0.  595 

-0.  507 

-0. 162 

0.  000 

3 

3.  5-5 

8 

0.028 

0.  208 

-0.099 

-0. 137 

0.000 

4 

5-8.  5 
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Category  II  Operations  (Cont. ) 


Ten  Minute  Forecast 

Order  (i) 

P1 

P2 

P3 

P4 

P. 

0 

Predictor 

Xi 

Limits  of 
Xi 

0 

0.  125 

0.  588 

0.  272 

0.  015 

0.000 

0 

1 

0.  G74 

-0.431 

-0. 230 

-0.013 

0.  000 

1 

1-2 

2 

-0.  125 

-0.  588 

-0.  251 

0.  126 

0.  839 

9 

40-80 

3 

-0.  125 

-0.  588 

-0.  255 

0.316 

0.  652 

8 

25-40 

4 

-0.  125 

-0. 578 

-0.  189 

0.  866 

0.  026 

7 

18.  5-25 

5 

-0.  125 

-0.  585 

-0.067 

0.  775 

0.  003 

(i 

12.  5-18.  5 

6 

-0.  115 

-0. 557 

0.  424 

0.  247 

0.  000 

f) 

8.  5-12. 5 

7 

-0.078 

-0.  219 

0.315 

0.054 

0.  000 

4 

5-8.  5 

8 

0.  194 

-0.  037 

-0.  145 

-0.  012 

0.000 

2 

2-3.  5 

Thirty  Minute  Forecast 

Order  (i) 

P1 

P2 

P3 

P4 

P5 

Predictor 

Xi 

Limits  of 
Xi 

0 

0.  184 

0.371 

0.375 

0.  066 

0.044 

0 

1 

-0. 184 

-0.  371 

-0. 297 

0.054 

0.  798 

9 

40-80 

2 

-0. 184 

-0.340 

-0. 278 

0.296 

0.  507 

8 

25-40 

3 

0.378 

-0.  101 

-0.  258 

-0.  024 

0.005 

1 

1-2 

4 

-0.  163 

-0.361 

-0. 215 

0.  650 

0.  089 

7 

18.  5-25 

5 

-0. 170 

-0.360 

-0.  131 

0.618 

0.043 

6 

12.  5-18.  5 

6 

-0.  125 

-0.275 

0.  118 

0.  255 

0.  027 

5 

8.  5-12.  5 

7 

0.  158 

0.  027 

-0.  166 

-0. 109 

-0.001 

2 

2-3.  5 

8 

-0.  053 

-0.  109 

0.013 

0.  140 

0.009 

4 

5-8.  5 

Sixty  Minute  Forecast 

Order  (i) 

P1 

P2 

P3 

P4 

P5 

Predictor 

Xi 

Limits  of 
Xi 

0 

0.  194 

0.  209 

0.  256 

0.  234 

0.  106 

0 

1 

-0.  168 

-0.  209 

-0.  136 

-0. 167 

0.  680 

9 

40-80 

2 

0.  250 

0.  066 

-0.  097 

-0. 147 

-0.072 

1 

1-2 

3 

-0.  135 

-0.  085 

-0.  146 

0.  110 

0.  256 

8 

25-40 

4 

-0.  194 

-0.  171 

-0.062 

0.489 

-0.  062 

6 

12.  5-18.  5 

5 

-0. 176 

-0. 181 

-0.  044 

0.425 

-0.024 

7 

18.  5-25 

6 

-0.  125 

0.008 

0.  071 

0.  104 

-0.058 

5 

8.  5-12.  5 

7 

-0.028 

-0.015 

0.  189 

-0. 102 

-0.044 

3 

3.  5-5 

8 

0.060 

0.094 

0.024 

-0.090 

-0.089 

2 

2-3.  5 
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Appendix  C 

Transition  Matrices 

The  following  one-step  transition  matrices  are  used  in  the  equivalent  Markov 
prediction  technique  and  can  be  interpreted  in  the  following  way: 


where 

=  P[SVRt  =  i  |  SVRo  =  j]  for  i,  j.  =  1.  2 . 5 

That  is,  P.j  is  the  transition  or  conditional  probability  that  SVR  at  time  t  will  be  in 
state  i  given  that  SVR  at  time  0  is  in  state  j. 
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1.  Category  I  Operations 


Two-minute  one-step  transition  matrix 


0.817 

0.152  0.030 

0.000 

0.000 

0.340 

0.426  0.219 

0.014 

0.000 

p  = 

0.025 

0.091  0.623 

0.260 

0.  000 

0. 004 

0.096  0.054 

0.  696 

0. 148 

_0. 000 

0.003  0.007 

0.  032 

0. 868_ 

Ten-minute  one 

-step  transition 

matrix 

0.  635 

0.275  0.092 

0.  000 

o. ooo- 

0.402 

0.318  0.234 

0.047 

0.000 

P  = 

0.034 

0.093  0.  594 

0.  277 

0.003 

0.011 

0.096  0.068 

0.  629 

0.  197 

0.000 

0.007  0.013 

0.  171 

0.  810 

2.  Category  II  Operations 

Two-minute  one-step  transition  matrix 


0.906 

0.083 

0.010 

0.  000 

0.  000~ 

0.178 

0.  624 

0. 198 

0.  000 

0.000 

P  = 

0.005 

0. 128 

0.  804 

0.063 

0.000 

0.  000 

0.000 

0. 104 

0.  891 

0.005 

0.  000 

0.000 

0.015 

0.  232 

0.753 

Five-minute  one-step  transition  matrix 

"5.  838 

0. 140 

0.  023 

0.  000 

o.  ooo” 

0.  207 

0.  574 

0.216 

0.  003 

0.000 

P  = 

0.011 

0. 138 

0.750 

0.  100 

0.000 

0.000 

0.002 

0. 130 

0.861 

0.  008 

0.  000 

0.  000 

0.018 

0.219 

0.  763_ 
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